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Abstract 

Where, how much and how efficiently the energy conversion takes place in a regenerator of a thermoacoustic engine are expressed 
using the axial distribution of acoustic work flow and heat flow. As a first step in determining the energy flows inside the regenerator, 
measuring methods of the work flow are briefly described and the experimental results in an acoustic resonator are shown. Applicability 
of these methods to the regenerator is discussed. 
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1. Introduction 

Coupled oscillations of pressure P, velocity U, and 
entropy s in an acoustic wave produce two energy flows; 
one is acoustic work flow I = A(PU) and the other is heat 
flow Q=T m S, where S = Ap m (sU) is hydrodynamic 
entropy flow, p m and T m are time averages of density and 
temperature of the gas, A is the cross-sectional area of a 
flow channel and, bars and angular brackets mean the time 
and the radial averages of the quantity inside of them [1^4], 
Furthermore, mutual conversions are possible between Q 
and /, when the acoustic wave travels in a small channel 
with an axial temperature gradient. Particularly when the 
channels are as narrow as those in a regenerator, the energy 
conversions are executed through a thermodynamic cycle 
similar to that in Stirling engines [5,6], As a result, it has 
become possible to eliminate pistons from Stirling engines 
by using the thermoacoustic technique [7]. 

Now many thermoacoustic devices working as piston¬ 
less Stirling prime movers and coolers [8,9] have been pro¬ 
posed as an extremely simple and reliable technology for 
the energy industry. On the other hand, the fact that acous¬ 
tic waves execute the energy conversions without pistons 
offers a new way to understand a heat engine, in which heat 
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flow Q and work flow I, instead of heat and work, constitute 
the basic concepts. In this paper, the energy conversions 
between Q and / are shown using energy flow diagrams 
in a thermoacoustic prime mover, and the experimental 
techniques for the measurement of / is introduced. 

2. Energy flows in thermoacoustic devices 

A thermoacoustic prime mover is schematically pre¬ 
sented in Fig. la from the conventional view point of a heat 
engine. A regenerator is the heart of the prime mover, since 
the energy conversion takes place inside of it. The heat 
exchangers at the sides of the regenerator act as heat baths 
with temperatures T H and T c . Input heat power Q H flows 
into the regenerator from the hot heat exchanger, while the 
heat power Q c flows out of the regenerator from the cold 
heat exchanger. The difference between Q H and Q c equals 
the output power A I. The efficiency of the energy conver¬ 
sion is given as A I/Qu and is always smaller than the Car¬ 
not efficiency (1 — T C /T H ), which is achieved only in an 
ideal regenerator without any irreversible processes. 

The axial distribution of the heat flow Q and the work 
flow I in the regenerator are shown in Fig. 1(b). Here the 
positive value of the energy flow represents the flow to 
the right, and the negative one represents to the left. The 
magnitude of the heat flow Q decreases (VQ< 0, where 
VQ means the divergence of Q) from Q H to Qc as it flows 
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Fig. 1. Schematic illustration of a Ihermoacoustic heat engine (aj and 
axial distributions of acoustic work flow and heat flow (b) and entropy 
flow (c). A gray horizontal line in (c) represents the entropy flow in an 
ideal regenerator, and a gray curve in (b) is the corresponding heat flow in 
the ideal regenerator. 


down the temperature gradient in the regenerator, while 
the work flow is amplified (V/ > 0) in the regenerator. 
From the first law of thermodynamics (energy conservation 
law), the relation 

V(0 + /)= 0 (1) 

always holds in the regenerator. 

The work source W _ V/ [3] represents the acoustic 
power production per unit length. In the conventional 
treatment of the heat engine, only the output power 
AI = Qh - Qc is discussed as a measure of the energy con¬ 
version, but from the thermoacoustic viewpoint, the mea¬ 
sure of the energy conversion can be locally addressed by 
W. In the case of / shown in Fig. 1(b), we see that the rate 
of the energy conversion increases with increasing .v. 

In addition to W, the increase of the entropy flow S per 
unit length, VS, is also important in characterizing the local 
energy conversion in the regenerator. The entropy flow S is 
given by the heat flow Q divided by the local mean temper¬ 
ature T m and is schematically shown in Fig. lc. In the 
steady state, VS is equal to the entropy production per unit 
length and unit time. The second law of thermodynamics 
states VS = 0 in the ideal regenerator, but otherwise 
VS > 0 in actual regenerators with irreversible processes. 
Hence, VS can be used as a measure of how far the thermo¬ 
dynamic processes associated with the local energy conver¬ 



Fig. 2. Acoustic resonator driven by an acoustic driver. The origin of the 
axial coordinate x is taken at the front of the driver. 


sion are away from reversible processes. In the case of S 
shown in Fig. 1(c), we see that the irreversible processes 
are more remarkable in the smaller x. Thus, we can see 
where, how much, and how efficiently the energy conver¬ 
sion takes place inside the regenerator from W and VS. 

To obtain an experimental understanding of thermoa¬ 
coustic energy conversions, it is essential to measure the 
axial distributions of Q and /, and determine W and VS. 
But the attempt to measure the energy flows has just begun. 
Measurement of the work flow / in a duct has recently been 
made possible, as will be introduced in the following sec¬ 
tion. The measurement of Q has not yet been established 
because of the difficulty of measuring the entropy oscilla¬ 
tion v. Tashiro has succeeded in the measurements of the 
temperature oscillation T with a thermocouple [11]. By 
using this technique, he plans to determine the total energy 
flow H= Q + /, since H is given as H = pC P (TU ) for an 
ideal gas, where C p denotes the isobaric specific heat of 
the gas. Acoustic heat flow Q can then be obtained as 
Q = H — I from the measured / and H. 

3. Work flow measurement 

3.1. Experimental methods 

Acoustic pressure P{x) = p(x)cxp( icot) with an angular 
frequency m is a function of the axial coordinate x, but it 
is independent of the radial coordinate r in the tube, when 
the tube radius is much smaller than the wavelength. The 
axial acoustic velocity U{x,r) = u(x,r)exp[i{cot + <p(x,r)}] 
depends on both x and r. Here the direction of the velocity 
is explicitly shown by its sign and we hereafter denote the 
velocity as U instead of U. By taking the radial average 
prior to the time average, the work flow l{x) = 
(p{x) U(x, r) ^ passing through the cross section at x is 
rewritten as 

I(x) = (P(x) V(x)}= jp(x)v(x) cos i//(x) (2) 

where V(x) = r(x)exp[i{<nf + i/flx)}] represents the radial 
average of the axial acoustic velocity. Therefore, the work 
flow / is derived from the measurements of P and V. 

Yazaki and Tominaga [10] demonstrated the acoustic 
power generation associated with spontaneous gas oscilla¬ 
tions from the measurements of the work flow I. They used 
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small pressure transducers for the measurement of the 
acoustic pressure, and a laser Doppler velocimeter (LDV) 
for the measurement of the axial velocity. This method 
was further applied to determine acoustic fields and acous¬ 
tic output powers in thermoacoustic engines [7,9,12-14]. 
The quality factor of an acoustic resonator was also deter¬ 
mined from the work flow measurements [15]. 

The measurement method of 7 which only involves 
pressure measurements has been reported by Fusco 
et al. [16]. In this method, P and V are calculated respec¬ 
tively from the sum and the difference of the pressure 
signals obtained by two adjacent pressure transducers. 
Simplicity of the measurements is a great advantage of 
this method, when compared to the method using the 
LDV. However, the estimated velocity and the work flow 
have not yet been verified yet by the simultaneous mea¬ 
surements of pressure and velocity. Next, we show how 
the work flow 7 can be determined in a simple acoustic 
resonator with these two methods in order to see the 
applicability of the latter method based on pressure mea¬ 
surements [17]. 

3.2. Work flow in an acoustic resonator 

The experimental setup is shown in Fig. 2. An acoustic 
resonator was made of a Pyrex glass tube with length 
L = 1.02 m and inner radius r 0 = 10.5 mm, and filled with 
atmospheric air. One end of the resonator is closed by a 
rigid plate, and the other was connected to a compression 
driver. 

The gas column in the resonator was driven at the fre¬ 
quency /= 298.5 Hz (the second mode), and the pressure 
amplitude at the closed end was kept at 1.00 kPa. The 
acoustic pressure P(x) = p(x)cxp{ ico t) was measured by a 
series of small pressure transducers mounted on the resona¬ 
tor wall via short ducts 10 mm in length. The axial acoustic 
velocity U c (x) on the central axis of the tube was measured 
using the LDV. Pressure and velocity signals were recorded 
with a multi-channel FFT analyzer. The measured velocity 
U c (x) — M C (x)exp [i{cot + tp c {x)}] was converted into the 
radial average V(x) = r(x)exp[i {cot + i/ffx)}] of the velocity 
using the theoretical solution of a laminar oscillating flow 
as 

v(x) = u c (x)/ 1.01, (3A) 

<A(x) = </> c (x) + 0-7- (3B) 

The phase angles of P{x) and U c (x) were accurately deter¬ 
mined by taking into account delays due to electrical cir¬ 
cuits and short ducts used in the pressure measurements. 

Fig. 3 shows the axial distribution of the work flow I(x) 
(open circles) derived from the simultaneous measurements 
of P and V. The measured I decreases from 7(0) = 24 mW 
to zero as it flows to the rigid end. The total input power 
7(0) from the compression driver gives the total energy loss 
E in the resonator. The experimental work source 
W= — HO)/L = — 24 W/m represents the dissipated power 



x (m) 

Fig. 3. Axial distribution of the work flow (open circles) derived from 
simultaneous measurements of P and V, (solid circles) from the pressure 
measurements with X = 0.2, and (solid triangles) with X = 0.6 m. 


per unit length in the resonator. The acoustic energy E s 
stored in the resonator was estimated as E s = 0.63 mJ from 
the axial distribution of P and V, and we obtained the qual¬ 
ity factor ci of the resonator as q = coEjE = 49 [15]. This is 
in good agreement with the theoretical value of 55 derived 
from a linear acoustic theory. 

3.3. Derivation of I from measured P 

The acoustic pressure P{x) and the radial average V(x) 
of the velocity at a position x can be estimated [16] from 
the sum and the difference of pressures P A = Fix — X/2) 
and Pn = P(x + X/2) measured with two pressure trans¬ 
ducers separated by a distance X. Namely, 


P( x ) - P a+F b 

W 2cos(LV/2) 

(4) 

V (x) — 1 (\ 

y> cop\ r 0 ) 2sin(LV/2) 

(5) 

where 


‘- t0 (' + 2 r, + 2 '>)’ 

(6) 


k 0 ( = m/a, a; adiabatic speed of sound) is the wave number 
in a free space, a and y are Prandtl number and the ratio of 
isobaric to isochoric specific heats of the gas, 5 means the 
boundary layer thickness given as S = yj2v/co with using 
the kinematic viscosity v. Eqs. (5) and (6) are valid when 
r 0 3><5 (boundary-layer approximation). The work flow 
I(x) is obtained by inserting P(x) and V(x) given by Eqs. 
(4) and (5) into (2). 

The work flow I(x) derived from the pressure data are 
shown as solid circles in Fig. 3. Here the distance X 
between the two measuring positions was chosen as 
X = 0.2 m. The relative error of the data obtained by pres¬ 
sure measurements was within 3%, and we conclude that 
this method is applicable to the resonator. 

The pressure measurements in this method are simpler 
and easier than the velocity measurements with the LDV. 
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Furthermore, there is no difficulty in the pressure measure¬ 
ments in stainless-steel tubes filled with pressurized gases, 
whereas the LDV requires transparent tubes for the trans¬ 
mission of the laser light. High-power thermoacoustic 
engines are usually composed of stainless-steel tubes 
because pressurized gases are used as a working gas. 
Hence, the measuring method of / only from the pressure 
signals is most suitable to the study of the acoustic powers 
in such engines [17]. 

It should be noted here that the choice of the distance X 
is important to accurately obtain I(x). We found significant 
error when X was taken as X = 0.6 m in the present resona¬ 
tor. The corresponding work flow I(x) is shown in Fig. 3 as 
solid triangles. The large error originates from the fact that 
the denominator in Eq. (4) tends to zero when X 
approaches the half wavelength. Therefore, the choice of 
X close to the half-wavelength should be avoided for 
the accurate evaluation of I(x) from the pressure 
measurements. 

3.4. Applicability to the work flow in a regenerator 

There are a few technical issues in determining of I from 
pressure signals in the regenerator. First, the validity of the 
boundary-layer approximation should be checked. Since 
viscous and thermal boundary layers of the gas fill up the 
channels in the regenerator, the boundary-layer approxi¬ 
mation is no longer applicable. Secondly, a very high reso¬ 
lution is needed in the phase measurements. The length of 
the regenerator is usually much smaller than the wave¬ 
length, and hence the phase difference between two pres¬ 
sure signals is small. On the other hand, the relative error 
of /becomes large when the phase difference is small. Thus, 
the derivation of / using only the pressure signals is attrac¬ 
tive, but challenging. 

Technical issues in measuring / in a regenerator with the 
LDV are that the regenerator should be made of transpar¬ 
ent materials and that its radius should be large enough 
compared to the size of the measuring volume of the 
LDV (typically, ^0.1 xO.l x 1 mm 3 ). In order to resolve 
these issues, we have constructed a liquid-piston Stirling 
engine [18] having a regenerator made of a glass tube. 
The working frequency of a liquid-piston Stirling engine 
is as low as 1 Hz, and the thermal boundary layer thickness 
becomes a few millimeters long for one-bar air. As a result, 
we could enlarge the channel radius of the regenerator to 
4 mm, which is much larger than those in other thermoa¬ 
coustic engines. We confirmed the applicability of the 
velocity measurement in such a narrow tube from prelimin¬ 
ary measurements [15]. Hence, we consider that the method 
with simultaneous measurements of P and V is more easily 
applied to the regenerator than that only with pressure 
measurements. The relatively large radius of the regenera¬ 
tor in the liquid-piston Stirling engine is also suitable for 
the temperature measurements using a thermocouple [11]. 


Measurements of P, V and T would be possible in the 
regenerator of the liquid-piston Stirling engine. 

4. Summary 

From the thermoacoustic viewpoint, the energy conver¬ 
sion inside the prime mover is expressed using the energy 
flows / and Q passing through the regenerator. Conse¬ 
quently, where, how much, and how efficiently the energy 
conversion takes place inside the regenerator are described 
by W=VI and VS = V(QlT m ). This is in clear contrast to 
the conventional viewpoint in which the quantities only on 
the outside of the regenerator is discussed. In order to dee¬ 
pen our understanding of heat engines from the thermoa¬ 
coustic viewpoint, measurements of / and Q in the 
regenerator is of great importance. 
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